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Drgania dwuwarstwowej lepkiej kropli 
Konebanua ABYXCIOUHOÓú BA3KOŃ KANJI 


The Oscillations of a Two-Layered Liquid Viscous Sphere 


§ 1. Drgania kropli cieczy idealnej i niescisliwej z uwzglednieniem 
sil grawitacyjnych rozwazal po raz pierwszy Lord Kelvin [1]. Podane 
przez niego rozwiązanie przeniósł Lamb [2] na przypadek drgań 
warstwy cieczy idealnej, pokrywającej sztywny rdzeń kulisty. W pracy 
[3] uogólnia się wspomniane wyniki na n warstw cieczy idealnych i nie- 
ściśliwych, pokrywających sztywny rdzeń, a w dyskusji uwzględnia sie 
przypadek kropli pokrytej warstwą innej cieczy. 

Drgania kropli cieczy lepkiej badał Chandrasekhar [4], a na- 
stępnie Reid [5] uogólnił jego wyniki na przypadek kropli ze sztyw- 
nym rdzeniem. Podali oni równania na częstości własne drgań tych 
kropli. | 

W naszej pracy wyprowadzamy .równanie na częstości własne drgań 
kropli cieczy lepkiej, pokrytej warstwą innej cieczy. Zagadnienie to 
obejmuje zarówno przypadek dyskutowany przez Chandrasekha- 
ra [4], jak również wyniki dla kropli dwuwarstwowej [3]. 

$ 2. Rozważamy, zgodnie z postawionym zagadnieniem, małe drgania 
ciekłej dwuwarstwowej kuli o promieniach powierzchni warstw równych 
odpowiednio R; i Rọ gdzie R< Ry. Przyjmujemy, że ciecz jest nie- 
ściśliwa, lepka i ciężka. Równania ruchu po odrzuceniu małych rzędu 
drugiego są następujące [4]. 

Ou — rę 
-pp = grad 0—v rot rot u (1) 
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div u=0 ; (2) 
=y +P (3) 
A8V=0, (4) 


gdzie u oznacza prędkość, v lepkość, o gęstość, 5p i 6V są zmianami 
ciśnienia i wewnętrznego potencjału grawitacyjnego, wywołanymi od- 
kształceniami kropli od kształtu kulistego. 

Równania odkształconych powierzchni r; = r, (9, p) i r» = r (0, q) 


piszemy w postaci szeregów funkcji kulistych y” (6,9) 
r" R, | toam r ©| (5) 


= Ry E tiam Y el, (6) 


gdzie RZ =at (t) dla «= 1,2 są współczynnikami zależnymi od czasu, 


a sumowanie rozciąga sie po m od —l do l i po l od 0 doo. Zmiana 
potencjału grawitacyjnego SV wyraża się następującymi wzorami [3]: 


dla To Srs 


a= 
= Dopey (gratin + Gree Ret of) YO. 


nor i ee, PE 
oraz dla OS r<rz (7) 
Pi, © P2— P a) > 
> GFT (z: =a Cm + Ro? a)” aaa 
Biorąc dywergencję memos (1) otrzymamy ze wzgledu na równanie (2) 
Aw = 0 (8) 
i ponieważ AV = 0, mamy: 
op 
+=0; (9) 
> p 
zatem: 
P- = Ż am APO FT dla r Lr (10) 
oraz: 


1+1 


R, 
en fet + pm me paja dla mZT=N (11) 
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gdzie Pp, P;, Pz są stałymi. Ze wzorów (3), (7), (10) i (11) otrzymujemy 
wyrażenia: 

ò 1 
o =—8V, + is = Że, |--o AYO ETAT rzesy e] yr (12) 
1 1,m 
dla nsr <r, 


3 
0 =V, H= = Ya Q+)at™ry”  dar<r (13) 
2 lm 


gdzie As, Ag, Az są stałymi. 
$ 3. Grad w jest wektorem biegunowym o funkcji skalnej [4]: 


1 
zaj =a (ar pi tt aim) >) dla r Sr <r 


aa © a AY q +1 dla r Ta; 


to znaczy, że zachodzą następujące równości: 
Ul + 1) 

(grad w), = "rst yp 

im Tr 
(grad o) = == 3 

m 

ży 1. oS dY, 

(grad a),= 2 Fan8 dr żę” 


> > 
uw i rot rot u są zatem także wektorami biegunowymi o funkcjach ska- 


(lym) 
larnych U (+ m (r) i odpowiednio — g ‘i WED gam, (14) 


u poprzez U (+) wyraża się w następujący sposób [4] 
u, =Je — opt ULT 1) A y em) y” 


„m 


(l,m) = 
A —G mt + 1 du ð £ 
145 2 AE 0 (15) 


Lm) gy” 
u, = >, e” %, mt 1 du‘ skie 
m 


rsin® dr op 
Po podstawieniu (12), (13), (14), (15) do (1), przyrównaniu współczynni- 
ków przy Y™ i opuszczeniu wskaźników otrzymamy: 


_ UT a +1 i 
> x U; wii U: > (a, T tA: > =| (16) 


104 Stefan Król, Maksymilian Piłat 


ý (1) 

o 
U; PU; = 
Ya Ya 


dU, +1) 
2 


dr Asr't? (17) 


Przyjęliśmy, że współczynniki a ® (t) w porównaniu (5), (6), (7), 10), (12) 
i (13) są postaci ade” 


Rozwiązaniami (16) i (17) są odpowiednio funkcje: 


1 1 (1) (1) 
1 ES a a. 1 
U(r) = Aor? Tia (qır) + Agr? aora] 07) Fa” Agr't? W LAŁ (18) 
ka a Desa í 
Valr)= Ayr? J, , 1 (qgar) +37 Ast , gdzie 
2 
>e D 6 
qı En , q2 Ye 
A a?) 

J (x) jest funkcją Bessela rzędu l. Stałe PJ = Aj, A, dla k = 2, 3, 4, 


5, 6, 7, wyznaczamy z nastepujacych wacabkóle granicznych: 


1. ur wyliczone z funkcji tworzacej U (m) i ze wzoru 


tR (1+ Xa 


2 yr) musza byé równe. 
im 


Stąd otrzymujemy dwa równania dla dwóch powierzchni: 
=T , u — r. , 
B= R, Š A r= Ro . 
2. Ciśnienie styczne P,e i p,, Musi znikać na powierzchni zewnętrznej. 
re Py, = O dla r = Ry. Prowadzi to do warunku: 


Lm a 1) U, 


u, 


+2 ey oT =0 dla r = Ry. 
3. Ciśnienie „Ło Pr, na powierzchni r = R; (1 +2 ajm y m) musi 


znikać. Ciśnienie to wyraża-się następującym wzorem: 


i Ur A 
P+ FP, + Spy + 2 rp ; gdzie 


2 a R 
pas mp?R,° (1) +4 (Pa — MIE IEA h -1) 


4. Prędkość radialna u, wyliczona z U;(r) powinna być równa w 
wyliczonej z U; (r) dla powierzchni r = Ry. di pzdr to do warunku 
U; = U; dla r = Ra. 
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5. Składowe prędkości prostopadłe do promienia powinny być ciągłe 


u ira ug |_ oraz u |: = ad dla r= R; 
dU; _ dU, AA 
Stąd ART dla r=R, 
6.. Napięcia styczne muszą być ciągłe w cieczy. 
FPr|- oraz 24 ls =p dla T= Ra 


stąd otrzymujemy dla r = Rz. następujące równanie: 


2 2 
pn {een _ 2 qd, +5 TN ee _ 2 dU> , ea] 


q? r dr dr? y, de dr? 


7. Warunek ciągłości ciśnień normalnych w cieczy prowadzi do rów- 
nania Pre| 4 = Prr| „ na powierzchni r = Rą (1 + Zaim YT”), 
które po przeksztalceniach przyjmuje postaé SR ve! MI 
(1) (2) 
U 
7 = ĝpa — 2Vap2 r 


Z wy aay warunków GSA następujący układ równań: 


Spy — 2V4P1 


Y ad, GB, = boat. 


k=1 
gdzie Ax są stałymi, o których była mowa poprzednio, natomiast 
A11 ~ Aza — As — Aig — Aga — a23 — Ang — Ag7 = Agg — Agi — Aga — agg — Ags — Aya — 
= 445 = A51 = Asg — Agi — a71 = a78 = 0, 


5 
A żę - 3 i 
as =R -J, +2 (4,R), a ie J_(142) (q,R,) 


1-2 
1 1 6 
8, =R, 1+1(,R,) oR Pe as TA ge 7, KUPI)” 
1 
6 Bee 


i 


3 
ag T [a+o+¿] R, TYE 1(q,R, )+3R, q Y 1 RIF 


MH 


+m) KU 1+1(,R,), 
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3 1 

ATE 1 
ayy = [+43] R 3 (4 49(0,8,) +R ŚŚ 2 dj 
RY 

= R2917"_ (14.2) (,R,), 


1+3 
a EA a,, = A e EATA 
25 == 
a,, =v,Ul+1) ES ; I +3) ORI —2R, * q,J’ SCJ 
ZE s 
ROZ J (14,2) (0,B) FAR, 29 (+ 04,8), 


E 2v (1+1) L E a (+2 
J 1 
eonia A EA 
1 
4 2 [2(0+1 R % 
Ass 3 Ri | 21+1 Pp, +, o) , a FRZ J1+1(q,R,), 
1 a RH 
=p? =— =e 
a,, = R; J- (+2) (9, R.), a Ry Ji+4(9,R,), as, — STH 
141 a 
234 raj e kabiny 7 ky m 2 y 
as 6 o RES 5 yA, J141 (,R,) +R, ag 14+7-@,R,), 
i 1 
r3 + 2 , 
as=2 P J-(41)0R)FR, 9,7, 1)(a, R.) 
4 1 
TES R : 
a,=— oF, Tipah J—R 97,1 E, 
FRL I „! BE 
>= Rz, as” — RE ye 


4 
an= pl (CENT sala? J44 1 (q,R,) + 3R, Ta, +2 (qR) + 


HR RD). 
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ko 


"z 


a, =P, Al (I+ 1+9/s] R, J 44,8, )+3R, J _ 44) GR) F 


R, q J”_ (+23) (a, R, 


M 


a, er raro +/s|e J +2 (q,R,) + 3R, » ay J’, +2 (9,R,) + 


M 


"2 r 
+R, % J LR), 


a 207D_py Rim _ 2041)" py rors 
5 A 6 2 


a 
a 
kJ 
= 


2 


21 1 2 
al A 81 


47 us 
21+1 R, (e, P,) ’ 


2,7 %910+D]3R, * HA - (q, R,)— 2R, q, ECU 


"w: 
R. *a,J’- (+2) (R), 


na 


a, =¥,p,10+1) se, -J-(+3) @R)—2 


=> a 
2 


a, =—¥,P,1(+1) [se, 


a,, 7-7, (L+ 1) |; 


J= (+7) (q, R) — 2R, tai aR), 


2v,1(1+1 = 
a |, 2 


dh 342%. BSE ay r + 27 (1+1)(1+2) 
O) e E a att Re a | 
4 R, 
T 2 
ass == FT pine Ph 


Po wyrugowaniu z tego układu równań stałych A; i= 1, 2, 


1 


sd 


otrzymujemy szukane równanie na częstości własne drgań dwuwarstwo- 


wej lepkiej kropli w postaci: 
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det | air | = 0, i, k=1, 2,..., 8. 


Znalezienie konkretnych rozwiązań tego równania wymaga stosowania 
metod numerycznych. 
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PESIOME 


B sro paore paccmMaTpUBAIOTCA TpaBUTaHNMOHHBIe MAJIbIe KOJIEÓAHNA 
ĄBYXCJIOŃHOŃ: KANJIN-- KOTODAA COCTOMT M3 MAKOTO Mapa, IOKPBITOTO 
CJIOeM MHOJi ZKMJĘKOCTM. IloJlaraeM, YTO STM ZKMĄKOCTM BASKU, TAK YTO 
rpoójremMa BKJIIOdaeT B cede Tak CJryuai KoJlreóaHni KANIA OAHopoamoú 
BA3KOŃ ZKMIKOCTM, KOTOpbIm paccmarpuBan HaHxHpacekxap [4], Kak M 
cryuaji KojreóaHnńi AByXCNOMHOM KANJI COCTABJIeHHOi M3 UACANBHBIX 
xnaKocreñ [3]. ; 

BbiBeqeHo ypaBHeHMe AJA YacToT COÓCTBEHHBIX KoJleÓaHuii. Ilony- 
ueHHoe ypaBHeHMe OKA3AJIOCb TPAHCLEHJEHTHBIM. 


SUMMARY 


The case of gravitational vibrations of a two-layered sphere is con- 
sidered. The vibrations system consists of a fluid sphere covered with 
a layer of another fluid. The fluids are supposed to be viscous so that 
our problem includes both, the vibrations of a homogeneous viscous 
sphere discussed by Chandrasekhar [4] and the vibrations of 
a two-layered sphere of ideal fluids [3]. 

The equation for the eigen-frequencies is obtained. This equation 
is transcendental. 
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